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Abstract 

The  importance  of  recognizing  that  each 
lower-order  aodel  used  for  rotorcraft  parameter 
identification  has  a  limited  range  of  applicabil¬ 
ity  Is  illustrated  in  some  detail.  Examples  are 
given  to  illustrate  the  use  of  conditioning  the 
test  input  signals  and  the  potential  of  using 
multi-axis  test  inputs  to  enhance  the  parameter 
ldentif lability.  The  paper  discusses  the  benefits 
and  limitations  of  using  frequency  sweeps  as 
flight-test  input  signals  for  identification  of 
frequency  response  for  rotorcraft  and  for  the 
subsequent  fitting  of  parametria  transfer-function 
models.  This  paper  demonstrates  the  major  role 
played  by  analytical  modeling  and  the  understand¬ 
ing  of  the  physios  involved  in  the  rotorcraft 
flight  dynamics,  particularly  understanding  the 
limit  of  lower-order  models,  in  achieving  success¬ 
ful  rotorcraft  parameter  identification. 


Notation  \ 

blade  coning  angle,  rad 

longitudinal  first-harmonic  flapping 
coefficient,  rad 

normal  acceleration,  ft/sec2 

lateral  cyclic  pitch,  rad 

lateral  first-harmonic  flapping  coeffi¬ 
cient,  rad 

longitudinal  cyclic  pitch,  rad 
stability  derivative  matrix 
control  derivative  matrix 
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rolling  moment  derivatives 
yawing  moment  derivatives 
roll  rate,  rad/sec  (or  deg/sec) 
pitch  rate,  rad/sec  (or  deg/sec) 
yaw  rate,  rad/sec  (or  deg/ sec) 

Laplace  transform  variable,  rad/sec 
time  constants,  sec 

longitudinal  airspeed  component,  ft/sec 
control  vector 

lateral  airspeed  component  (also  inflow 
velocity),  ft/sec 

airspeed,  ft/sec 

vertical  airspeed  component,  ft/sec 

state  vector 

side  force  derivatives 

quasi-static,  vertical  damping,  1/sec 

quasi-static  oolleetive  effectiveness,  _ _ 

ft/see2/rad  /" 

*4- 

pitch  attitude,  rad  S 

roll  attitude,  rad 
tail-rotor  collective  pitch  rad 
collective  pitch,  rad 

lateral  stick  displacement  (also,  1 ' 

‘AS>>  !L_ 

—  aft,  a 


\fi'l 


□□ 


&rp  rudder  pedals  displacement,  in. 

tQ  collective  stick  displacement,  in. 

«p  rudder  aotuator  displaceaent,  deg 

6,  aileron  actuator  displaceaent ,  deg 

t  tine  delay,  sec 

a  standard  deviation 

B  side  slip  (also  flapping  angle),  rad  (or 

deg) 

undamped  natural  frequency,  rad/sec 
C  damping  ratio 

eigenvalues 


Introduction 

For  many  years,  research  has  been  conducted 
by  the  helicopter  industry  and  government  agencies 
to  develop  global  computer  models  such  as  C-8l , 1 
REXOR , 2  and  CAMRAO,^  which  include  comprehensive 
aerodynamics,  structural  dynamics,  engines,  and 
flight-control  systems.  For  flight  dynamics  and 
control  applications,  the  validity  of  these  inter¬ 
disciplinary,  global  computer  codes  can  only  be 
established  by  comparing  their  predicted  responses 
with  those  actually  measured  in  flight.  Because 
of  the  numerous  assumptions  and  approximations 
used  and  the  large  number  of  parameters  involved 
in  a  global  computer  code,  it  is  extremely  diffi¬ 
cult  to  fully  establish  its  credibility  by  adjust¬ 
ing  these  assumptions  and  parameters  to  fit  the 
measured  data. 

Simplified,  special -purpose  analytical  models 
that  are  easier  to  operate  and  to  comprehend  can 
satisfy  the  needs  in  dealing  with  many  aspects  of 
flight-dynamic  problems  each  as  helicopter  flying- 
qualities  evaluation,  design  of  stability  and 
control  augmentation  systems  (SCdS),  and  ground 
simulator  validation.  These  siaplificd  analytical 
models  nay  be  directly  developed  from  the  flight- 
test  data  by  using  modem  system  or  parameter- 
identification  techniques. 

Significant  advancement  has  been  made  in  the 
field  of  aircraft  state  estimation  and  parameter 
identification  in  the  last  two  decades.  Two  sym¬ 
posia and  an  MUM)  lecture  series6  have 
reviewed  the  state  of  technology  in  this  field. 
These  symposia  and  lecture  series  indicated  that, 
while  significant  progress  has  been  made  in  the 
Identification  of  the  stability  and  control 


parameters  of  the  fixed-wing  alrcrat,  by  compari¬ 
son,  progress  has  been  relatively  slow  in  its 
rotary-wing  counterpart.  Several  important  issues 
and  problems  are  still  facing  rotorcraft  parameter 
identification  which  need  to  be  resolved  if  major 
advancements  are  to  be  achieved. 

Unlike  the  flight  dynamics  of  fixed-wing 
airoraft,  rotary-wing  aircraft's  are  characteris¬ 
tically  thosa  of  a  high-order  system.  Significant 
intar-axla  couplings  axist  fcr  single  main-rotor 
helicopters^;  dynamic  interactions  are  present 
between  the  engine  and  drive  train/rotor  system8 : 
and  high-order  effects  such  as  rotor  dynamics^-12 
and  Inflow  dynamics^-18  are  inherently  present  in 
the  system.  The  large  number  of  degrees  of  free¬ 
dom  associated  with  the  coupled  high-order  dynam¬ 
ics  leads  to  a  large  number  of  unknown  parameters 
that  have  to  be  identified,  making  it  extremely 
difficult  to  achieve  a  successful  application  of 
system  or  parameter-identification  techniques. 

Special  purpose  lower-order  models,  valid  for 
a  limited  range  of  frequencies  or  input-output 
magnitudes,  may  be  used  to  reduoe  the  number  of 
unknown  parameters  to  be  identified.  The  selec¬ 
tion  of  a  lower  order  model  is  often  dictated  by 
its  usage.  For  flight-control  applications,  for 
exmaple,  a  model  valid  for  a  higher  frequency 
range  la  required  for  the  design  and  analysis  of 
high  bandwidth  flight-control  systems like¬ 
wise,  a  model  valid  only  for  a  lower  frequency 
range  such  as  one  of  quasi-static  rigid-body 
dynamic  models  can  be  quite  adequate  for  low  band¬ 
width  flight-control  systems.  Given  a  frequency 
range  of  interest,  the  determination  of  an  appro¬ 
priate  dimensionality  (or  system  order)  for  a 
lower-order  helicopter  model  is  not  a  simple  task, 
however.  Generally,  in  mathematics,  the  higher 
the  order  of  a  model,  the  better  the  fit  to  the 
observed  input-output  data  over  the  range  of 
frequencies  of  interest.  What  role,  then,  is  to 
be  played  by  physics  and  what  by  mathematics  (or 
identification  algorithms)?  What  type  of  a  lower- 
order  model  ehould  one  use,  parametric  or 
nonparametrio? 

Whether  a  parametric  or  a  nonparametrio  model 
Is  to  be  used  for  identification  depends  on  the 
usage  of  the  model  and  the  identification  proce¬ 
dure  selected.  For  flight  dynamics  and  control 
applications,  parametric  models  such  as  state- 
space  models  or  transfer-function  models  are  pre¬ 
ferred  over  nonparametrio  models  such  as  impulse 
or  frequency  response.  Identification  of  a  param- 
etrio  model  can  be  conducted  in  either  the  time 
domain1’’2^  or  the  frequency  domain.28  The  latter 
can  also  be  employed  for  identification  of  a  non- 
parametrie  modal  first,  then  a  parametric 
model. 2“-30  sophisticated  software  packages  for 
both  time-  and  frequency-domain  system  (or 


parameter)  identification  have  become  n a— arc tally 
available  in  the  past  few  years.31,32  The  Icey  to 
a  successful  rotor craft  parameter  identification 
My  lie  in  the  selection  of  appropriate  loner- 
order  aodels  discussed  above ,  and  in  the  design  of 
proper  test  inputs."*35 

The  design  of  control  input  signals  for 
flight  testing  is  crucial  for  the  identification 
of  loner-order  helicopter  flight-dynaaics 
Models.  Test  inputs  nay  be  designed  to  optlaise 
the  infomation  content  in  the  helicopter  response 
aeasurenents  within  the  frequency  range  of  inter¬ 
est  and,  at  the  sane  tine,  nininize  the  ezeitation 
of  the  helicopter  nodes  outside  the  frequency 
range  of  the  nodel  valid  for  its  intended  usage. 
Helioop ter  responses  to  a  control  input  can  be 
considerably  nonlinear,  especially  in  the  vertical 
axis. 15  Care,  therefore,  nust  be  ezeroised  to 
devise  the  size  of  the  test  inputs  so  that  the 
Magnitude  of  the  helicopter  response  Is  eomnensu- 
rate  with  the  intended  use  of  the  identified 
Model.  TO  naxinize  the  inforMtlon  content  in  the 
helicopter  response  data  and  thereby  enhance  the 
identif  lability  of  the  unknown  parameters ,  espe¬ 
cially  those  associated  with  internals  coupling, 
it  nay  be  desirable  to  eaploy  Multi-axis,  instead 
of  single-axis,  test  inputs. 

Other  difficulties  associated  with  rotorcraft 
parameter- Identification  include  the  relatively 
high  noise  level  contamination  in  the  response 
Measurements,  and  the  inherent  instability  of  the 
basio  aircraft.  Means  of  overooalng  seas  of  these 
difficulties  will  be  discussed  in  conjunction  with 
the  design  of  test  inputs  in  the  remainder  of  the 
P*P*r. 


Modeling  lelioopter  Dynamics  for 

It  is  important  to  recognize  that  each  ana¬ 
lytical  Model  used  to  fit  the  finite  input-output 
data  has  a  Halted  range  of  applicability.  Recog¬ 
nizing  this  fact  is  of  special  laportanee  in  Mod¬ 
eling  helicopter  dynaeios  for  parameter  identifi¬ 
cation,  because  a  he lloop ter  has  a  relatively 
large  number  of  degrees  of  freedom.  Awareness  of 
the  llaitations  of  the  analytical  aodel  is  also  of 
vital  laportanee  in  the  design  of  flight-test 
experlaents,  in  processing  the  test  data  and 
selecting  of  parameter  identification  algorithms, 
and  in  correlation  and  interpretation  of  the  iden¬ 
tification  results. 

The  selection  of  the  range  or  applicability 
of  the  nodel  is,  of  course,  dictated  by  the  Moti¬ 
vation  behind  the  identification.  For  example,  if 
the  motel  is  to  be  used  to  validate  a  ground 
siaulation  involving  large  Maneuvers,  then  the 


appropriate  candidate  nodel  would  be  nonlinear, 
not  a  linear  nodel.  If,  on  the  other  hand,  the 
motel  is  to  be  used  for  correlation  with  flying- 
qualities  evaluation  results  or  for  fllght-oontrol 
system  design  involving  possible  use  of  high 
gains,  then  a  better  candidate  would  be  a  linear 
time-invariant  nodel,  rot  a  nonlinear  aodel.  Two 
lower-order,  linear,  time-invariant  models  are  now 
exanined  to  Illustrate  their  validity  in  terms  of 
frequency  range. 

Validity  and  Limitation  of  a  Quasi-Static 
Slx-Dearee-of -Freedom  Rigid-Body  Model 


Consider  a  linear,  time-invariant,  nine  DOF 
representation  of  the  helicopter  dynamics,  which 
includes  six  DOF  rigid-body  dynamics  and  three  DOF 
tip-path  plane  dynamics  for  flapping  motion.  This 
node!  is  useful  for  the  design  of  high  bandwidth 
flight-control  systems.10-11  The  nodel  has  the 
following  form: 
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representing  the  state  vector  of  the  rotor  and  the 
body  dynamics;  u  «  (A,  B,c  «c  #,.)•  is  the 
rotor  oontrol  vector;  Fi,  and  Gl  are  matrices  of 
stability  and  oontrol  derivatives.  Table  1(a) 
Chews  a  set  of  numerical  values  of  the  stability 
and  oontrol  derivatives  of  a  representative  tee¬ 
tering  rotor  helicopter  in  hover.  The  representa¬ 
tion  is  a  lOth-order  system  (the  decoupled  coning, 
yaw,  and  vertical  degrees  of  freedom  are  not 
Included) . 


A  simplified  six  DOF  model  which  assumes 
instantaneous  response  of  the  rotor  nay  be  devel¬ 
oped  from  the  nine  DOF  motel  (1)  using  the  steady- 
state  equation  for  the  rotor  dynamics .  The  resul¬ 
tant  quasi-static  six  DOF  rigid  body  dynamics 
motel  has  the  following  form: 

*T  “  <F22  ‘  F21F11F12)kF  *  (C2  "  F21P11<VB 

The  reduced  lower-order  nodel  (6th-order)  is  shown 
in  Table  1(b).  The  two  low-frequency  complex 
modes  are  the  same  for  both  models,  but  the  real 
roots  are  more  widely  separated  for  the  lower- 
order  motel.  The  differences  between  these  two 
motels  are  more  olearly  shown  in  their  frequency 


A  comparison  of  the  frequency  response  of  the 
fuselage  pltoh-and-roll  rate  to  the  lateral  cyclic 
input  for  the  10th -order  system  and  the  6th -order 


system  as  given  in  Tables  1(a)  and  1(b)  are  shown 
in  Figs.  1  and  2.  It  is  worth  noting  from  these 
figures  that,  although  the  character 1st io  frequen¬ 
cies  of  the  rotor  dynamics  and  the  body  dynamics 
are  well  separated  (by  sore  than  a  factor  of  10), 
the  frequency  range  of  applicability  of  the  sim¬ 
plified  6th-order  model  is  good  only  up  to  approx¬ 
imately  twice  that  of  tbe  highest  rigid-body  char¬ 
acteristic  frequency.  In  fact,  if  the  amplitude 
ratio  miss-match  of  a  6  dB  or  phase  differential 
of  12*  (whichever  is  more  stringent),  is  used  as 
the  criterion,  the  maximum  frequency  that  satis¬ 
fies  this  criterion  is  approximately  2.5  r/s. 

Thus,  it  cannot  be  overemphasized  that  a  simpli¬ 
fied  six  DOF  analytical  model  such  as  Eq.  (2)  has 
a  limited  frequency  range  of  applicability. 

To  further  demonstrate  the  point,  consider 
two  longitudinal  cyclic  inputs,  one  "doublet"  with 
the  dominant  frequency  content  of  approximately 
1.57  r/s  which  is  within  the  frequency  range  of 
applicability  of  the  6th-order  model,  and  the 
other  input  with  a  dominant  frequency  content  of 
9  r/s  which  is  beyond  the  frequency  range  of  its 
applicability.  Figures  3  and  4  show  a  comparison 
of  the  6th-  and  lOth-order  models  using  the  two 
inputs  mentioned  above.  These  figures  clearly 
indicate  the  importance  of  recognizing  the  fre¬ 
quency  range  of  applicability  of  a  simplified 
analytical  model  in  designing  and  conditioning 
test  inputs.  Input  design  will  be  disoussed  later 
in  the  paper. 


Models  in  Hover 


Mow  consider  a  linear  time- invariant  three 
DOF  representation  of  the  helicopter  dynamics, 
which  includes  dynamic  inflow,  flapping,  and  the 
vertical  motion  of  the  helicopter1**  in  hover. 
Again,  this  model  may  be  useful  for  flight-control 
system  design  involving  possible  use  of  high  gains 
in  the  vertical  axis.  A  set  of  stability  and 
control  derivatives  of  an  articulated  rotor  heli¬ 
copter1”  for  this  4th -order  system  is  shown  in 
Table  2.  Following  the  procedure  used  to  reduce 
Eq.  (1)  to  Eq.  (2),  two  reduced-order  analytical 
models,  one  without  inflow  dynamics,  and  the  other 
without  inflow  and  flapping  dynamlos  are  calcu¬ 
lated,  with  the  results  shown  also  in  Table  2. 
These  two  reduced  order  models  are  represented  by 
a  3rd-order  and  a  Ist-order  system,  respectively, 
the  latter  being  the  conventional  (quasi-static) 
vertical  dynamics, 

■  3  V  +  \*o  (3> 

with  Zy  «  -0.304  (1/sec),  and 

2.  <  -293.34  ft/a2/rad. 
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To  provide  some  background  for  the  results  to 
be  discussed  in  the  subsequent  sections,  it  is 
important  to  compare  the  calculated  vertical- 
acceleration  responses  from  the  three  models  in 
Table  2  with  that  measured  from  flight  using  a 
CH-47B  research  aircraft. 3”  Shown  in  Fig.  5  is  a 
typical  vertical  acceleration  response  to  collec¬ 
tive  step  input  with  a  magnitude  of  0.62  in.  (or 
0.0201  rad  collective  pitch  change).  The  data 
were  taken  initially  with  a  sampling  frequency  of 
107.75  Hz,  then  passed  through  a  5  Hz  filter  to 
remove  the  3/rev  (11  Hz)  and  higher  harmonic 
vibratory  noise.  The  filtered  data  were  then 
decimated  by  a  factor  of  5  to  yield  21.55  Hz  data, 
as  plotted  in  Fig.  5.  As  can  be  seen,  the  tran¬ 
sient  is  characterized  by  an  oyershoot,  which 
immediately  follows  the  abrupt  change  in  collec¬ 
tive  pitch.  The  response  of  the  rotor  RPM  was  not 
recorded,  but  some  small  variations  in  rotor  RPM 
on  the  order  of  35,  were  observed  from  a  video 
recording  of  the  oookpit  instrumentation  during 
the  flight.  The  response  of  the  RPM  is  slower 
than  that  of  normal  acceleration,  typically  droop¬ 
ing  down  to  its  minimum  at  about  1  sec  after  an  up 
collective,  and  the  transient  lasts  for  about 
3  see  before  reaching  steady  state. 

The  calculated  perturbation  of  the  vertical 
acceleration  responses  to  a  0.62  in.  step  input  in 
collective  pitoh  is  shown  in  Fig.  6  for  the  three 
models  shown  in  Table  2.  The  4th -order  model, 
which  Includes  dynaslc  inflow,  flapping,  and  the 
vertical  motion  of  the  helicopter,  matches  the 
flight  data  reasonably  well.  The  3rd  order  model, 
which  neglects  the  inflow  dynamics,  results  in  a 
degraded  natch  with  the  flight  data  in  the  initial 
transient;  and  the  quasl-statie  Ist-order  model 
completely  falls  to  capture  the  initial  overshoot, 
which  is  a  result  of  the  combined  effects  of  the 
two  higher  frequency  modes,  the  inflow  mode  and 
the  flapping  mode.18 

Tbe  validity,  in  terms  of  frequency  range,  of 
the  two  lower-order  models  can  be  seen  more 
clearly  by  comparing  their  frequency  responses 
with  those  of  the  4th-order  model  as  shown  in 
Fig.  7.  The  calculations  cover  the  frequency 
range  of  0.1  to  100  rad/sec.  A  resonant  peak  at 
the  frequency  of  about  17  rad/sec  is  evident  from 
the  4th -order  model.  This  is  attributable  to  the 
destabilized  flapping  node  induced  by  the  dynamic 
inflow18  and  can  clearly  be  seen  by  comparing 
Fig.  7  (a)  and  (b).  The  frequency  range  of  appli¬ 
cability  of  the  Ist-order,  quasi-static  model  is 
considerably  lower  compared  to  the  3rd -order 
model.  If  the  criterion  described  previously  for 
the  six  DOF  rigid-body  model  is  applied  here,  the 
frequency  range  of  applicability  is  less  than 
5  r/s.  The  transfer  functions  of  the  vertloal 
acceleration  to  the  oolleotive  stick  displacement 
are  shown  in  Table  3.  It  is  Interesting  to  note 


that  both  tho  4th-order  and  the  3rd-order  Models 
exhibit  nonalnlnun  phase  characteristics. 


Selection  of  Teat  Incuts 

The  test  Input  Is  one  of  the  aost  Important 
factors  affecting  the  accuracy  with  which  Model 
paraaetera  can  be  identified .  In  fact ,  for  a 
specified  aeasureaent  system  there  is  a  theoreti¬ 
cally  achievable  aaxiauM  parameter  identification 
accuracy  corresponding  to  each  flight  test  input 
used  to  exoite  the  helicopter.  In  the  design  of 
flight-test  inputs,  two  important  questions  oust 
be  answered:  (1)  What  should  the  input  function 
be  so  that  all  the  helicopter  dynaaie  aodes  within 
the  frequency  range  of  interest  are  properly 
excited?  (2)  For  how  long  should  the  data  record 
be  taken  to  enable  identification  of  the  paraa- 
eters  to  a  desired  level  of  accuracy? 

Many  tiae-doaain  and  frequency-dcaain 
approaches  to  the  design  of  flight-test  inputs  for 
aircraft  parameter  identification  have  been  exten¬ 
sively  discussed  in  the  literature. These 

methods  invariably  involve  maximisation  of  seme 
function  of  Fisher's  information  matrix  or  minimi¬ 
sation  of  some  function  of  the  oovar lance  matrix 
of  the  parameter- identification  error.  Sinoe  all 
the  input-design  procedures  begin  with  a  set  of 
a  priori  values  of  the  very  parameters  that  are  to 
be  identified,  no  input-design  procedure  can  truly 
be  optimal  in  practice.  For  time-domain  identifi¬ 
cation  of  aircraft  stability  and  control  param¬ 
eters,  some  practical  approaches  to  the  design  of 
"multi-step"  test  Inputs  have  been  developed  using 
frequency  analysis^  or  Walsh  functions. 3* 

Figure  8  shows  same  typical  multi-step  input 
signals  and  their  frequency  distribution  plots. 

For  oamparison  purposes,  a  pulse  input  and  a 
smooth  test  signal  designed  using  the  method  of 
Ref.  33  are  also  shown  in  the  figure.  For  the 
pulse  or  the  multi-step  input  signals,  there  is 
considerably  more  power  in  the  high-frequency 
region  because  of  the  sharp  edges  associated  with 
these  signals.  Proper  conditioning  of  these 
signals  is,  therefore,  required  if  a  successful 
identification  of  lower-order  models  is  to  be 
achieved. 

Conditioning  of  Innut  Signals 

Because  of  the  limited  frequency  range  of 
applicability  for  lower-order  models,  proper 
design  of  flight-test  inputs  for  identification  of 
stability -and -control  parameters  of  the  helioop ter 
is  even  more  important.  The  highest  frequency 
content  of  the  flight-test  inputs  should  be  prop¬ 
erly  controlled  so  as  to  be  consistent  with  the 
frequency  range  of  applicability  of  the  model. 
Otherwise,  the  results  of  parameter  identification 


would  be  very  difficult  to  interpret  and  to  corre¬ 
late  with  those  generated  from  a  more  complete 
model.  Experience  in  applying  the  method  outlined 
in  Ref.  34  to  a  variety  of  aircraft  has  shown  that 
the  dominant  frequencies  contained  in  the  test 
input  signals  are  all  lower  than  the  highest  char¬ 
acteristic  frequency  of  the  model,  thus  assuring 
the  consistency  of  the  input  frequency  distribu¬ 
tion  and  the  frequency  range  of  applicability  of 
the  model.  However,  to  insure  that  the  high  fre¬ 
quency  content  of  these  input  signals  does  not 
exceed  the  frequency  range  of  applicability  of  the 
model  selected  for  parameter  identification,  a 
low-pass  filter  may  be  used  to  operate  on  the 
pilot  input  or  on  an  automatic  input  device  as 
schematically  shown  in  Fig.  9. 

The  validity  of  this  concept  may  be  seen  by 
comparing  Fig.  10  with  Fig.  4  as  discussed  in  a 
preceding  section.  In  Fig.  10,  the  high-frequency 
(9  r/s)  doublet  input  applied  in  Fig.  4  was  prop¬ 
erly  shaped  using  a  2nd -order  filter  with 

»  2.5  r/s  and  c  =  0.7.  To  achieve  a  better 
oamparison,  the  amplitude  of  the  doublet  in  Fig.  4 
is  increased  so  that  the  filtered  input  in  Fig.  10 
has  the  same  peak  amplitude  as  that  in  Fig.  4. 
Because  of  the  attenuation  of  the  high-frequency 
content  of  the  input,  a  significant  improvement  is 
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Multi-axis  Test  Inputs 

Multi-axis  test  inputs  nay  be  designed  to 
enhance  the  identif lability  of  the  helicopter 
unknown  parameters,  particularly  those  associated 
with  the  interaxis  coupling.  For  example,  the 
identification  of  the  lateral-directional  stabil¬ 
ity  derivatives  of  a  large  fixed-wing  cargo  air¬ 
plane  nay  be  considerably  improved  with  test 
inputs  involving  simultaneous  use  of  aileron  and 
rudder.  Shown  in  Table  4  are  the  separate  aileron 
inputs  and  rudder  inputs,  both  with  a  doublet  and 
a  pulse  having  their  dominant  frequency  tuned  to 
that  of  the  Dutch-roll  node.  For  this  particular 
example,  the  aileron  Inputs  are  much  more  effec¬ 
tive  than  the  rudder  inputs  in  achieving  smaller 
errors  in  identification  of  the  lateral- 
dlreotional  stability  derivatives.  The  alroraft 
responses  to  these  inputs  are  not  exoessive  and 
are  on  the  sane  order  of  magnitude.  Figures  11 
and  12  show  the  aircraft  responses  to  the  rudder 
doublet  and  the  aileron  doublet,  respectively. 

The  practical  Walsh-funetlon  method  of 
Ref.  34  was  then  used  to  design  test  inputs  with 
simultaneous  use  of  aileron  and  rudder.  Two  sets 
of  "suboptimal"  inputs  are  shown  in  Table  5.  It 
is  evident  from  Tables  4  and  5  that  for  the  sane 
data  length  of  16  seo,  the  suboptiaal  inputs  pro¬ 
vide  better  identification  results  than  the  con¬ 
ventional  pulse  and  doublet  inputs.  Because  of 


the  simultaneous  use  of  both  tho  aileron  and  the 
the  subopt imal  Inputs  also  permit  simulta¬ 
neous  identification  of  all  the  control  deriva¬ 
tives  in  addition  to  all  the  stability  deriva¬ 
tives.  Thus,  they  are  effective  and  economical. 
Figure  13  shows  the  aircraft  responses  to  the 
suboptlsal  Input  II  of  Table  5. 

For  oospar Ison  purposes,  an  evaluation  was 
also  node  for  an  alternate  set  of  rudder  and 
aileron  inputs  devised  years  ago  by  flight 
research  engineers.  This  set  of  Inputs  can  easily 
be  performed  by  a  test  pilot  by  slaply  observing 
the  bank  angle.  As  shown  in  Table  6,  a  rudder 
step  is  first  applied.  When  the  bank  angle 
reaches  a  prescribed  value,  say  20*,  an  aileron 
step  is  then  applied  to  return  the  bank  angle  to 
0*.  At  this  point  the  rudder  step  is  taken  out, 
leaving  in  the  aileron  input.  The  aileron  is 
later  reversed  to  bring  the  bank  angle  back  to 
0*.  As  shown  in  Table  6,  this  set  of  praotieally 
iwplown table  inputs  by  the  pilot,  though  not  as 
effective  as  the  subopt inal  inputs  of  Table  5,  is 
significantly  better  than  either  the  conventional 
aileron  inputs  or  rudder  inputs.  Figure  14  shows 
the  simulated  aircraft  response  to  this  set  of 
praotical  lateral -directional  inputs. 

Another  practical  input  signal  that  has  been 
used  successfully  in  a  nuwber  of  rotoreraft  flight 
fcesfcd®?»3®»39,40  nonrotorcraft  flight 
tests” * 3  is  the  "frequency  sweep."  The  fre¬ 
quency  sweep  provides  good  oontrol  over  the  fre¬ 
quency  range  of  interest  and  so  lends  itself  natu¬ 
rally  to  frequency-domain  identification. 


TWo  typical  concatenated  lateral  stick  (A^j) 
frequency  sweeps  completed  during  the  hover  tests 
of  a  teetering  rotor  helicopter^?  are  shown  in 
Fig.  IS.  These  tests  used  pilot-generated  rather 
than  computer-generated  Inputs.  The  sweep  is 
initiated  with  two  low-frequency  sinusoidal  shaped 
cycles,  with  the  periods  corresponding  to  the  low- 
frequency  bound  of  the  desired  identification 
range.  For  these  particular  tests,  a  nominal 
identification  range  of  0.4-15.0  rad/sec  was 
desired,  so  the  low-frequency  period  was 
T  ■  16  see.  After  the  initial  two  low-frequency 
oyoles,  the  oontrol  is  moved  at  gradually  increas¬ 
ing  frequency  for  another  50  seo.  By  the  end  of 
the  run,  the  oontrol  is  being  driven  quite  rapidly 
(about  4  Ha  in  Fig.  15)  with  generally  smaller 
displacements.  The  oontrol  is  then  returned  to 
trim,  ending  the  approximately  90  sec  test 
period.  At  least  two  repeat  sweeps  are  performed 
oonseoutlvely  for  each  oontrol  to  ensure  a  suffi¬ 
cient  amount  of  dynamic  data  for  the  identifica¬ 
tion  prooess. 


The  input  autospectrum  (Gs<(w))  in  Fig.  16 
displays  the  frequency  distribution  of  the  concat¬ 
enated  excitation  signal  (A^  of  Fig.  15).  The 
frequency-sweep  produces  nearly  oonstant  power  in 
the  range  of  0. 3-7.0  rad/ sec.  The  spectral  con¬ 
tent  apparent  below  the  minimum  average  input 
frequency  of  w  *  0.4  rad/sec  (T  *  16  seo) 
results  from  the  various  nonsinusoidal  low- 
frequency  input  signal  details.  At  high  frequenoy 
<«  >  7  rad/sec),  the  reduced  autospeotrum  refleets 
the  deliberate  reduction  in  input  amplitude  (this 
is  also  apparent  in  Fig.  15)  which  was  neoessary 
to  avoid  overly  exciting  the  rotor-pylon  dynam¬ 
ics.  The  pilots  could  comfortably  generate 
sizable  inputs  up  to  a  frequency  of  about  4  He. 

The  frequency-sweep  input  is  especially  well 
suited  for  frequency-domain  identification  proce¬ 
dures  because: 

(1)  The  input  autospeotrum  is  generally 
constant  over  the  desired  frequenoy  range. 

(2)  The  wave  form  is  roughly  symmetric,  so 
the  mean  and  linear  drift  oan  be  extracted  by 
noniterative  algebraic  means— a  requirement  for 
nonpar ana trio  identification. 

(3)  The  input  and  output  wave  forms  are 
smooth  and  regular,  so  the  resulting  apeotral 
functions  are  well  behaved.  (Miltl-step  inputs, 
if  not  conditioned  properly,  are  not  suitable  for 
frequenoy -domain  identification  because  the  abrupt 
wave  form  generates  undesirable  side-lobes  in  the 
spectral  functions.) 

(4)  The  frequency  buildup  is  slow  and  steady 
(as  seen  in  Fig.  15)  which  allows  the  sweep  to  be 
terminated  if  structural  node  excitation  is 
enoountered.  Multi-step  inputs  may  "ring"  the 
aircraft  natural  modes  simultaneously  unless  con¬ 
ditioning  is  used. 

Some  problems  encountered  with  frequency- 
sweep  testing  are: 

(1)  The  long  testing  run  needed  to  identify 
low-frequency  dynamics  generally  requires  that  the 
hover  tests  be  conducted  with  the  stability  aug¬ 
mentation  system  engaged.  Thus,  the  extraction  of 
open-loop  dynamics  is  sensitive  to  correlation  of 
the  exoltation  signal  (e.g.,  surface  deflection) 
with  the  feedback  signal.  Success  in  extracting 
open-loop  dynamics  from  closed-loop  testing  is 
essentially  a  signal-to-noise  problem  that  oan  be 
overcame  if  the  feedback  gains  are  not  very 
high. 

(2)  For  highly  ooupled  (e.g.,  single-rotor) 
helioopters,  off-axis  excitation  is  unavoidable 
and  must  be  regulated  by  some  moveewnt  of  the 
"nonswept"  controls.  (This  diffloulty  is 


particularly  noticeable  in  the  hover  oondition.) 
Such  secondary  regulation  is  acceptable  and  oan  be 
accounted  for  in  the  analysis  if  the  swept  and 
nonswept  inputs  are  not  fully  correlated. 

(3)  As  seen  in  Fig.  15,  the  pilot's  input 
amplitude  changes  as  a  function  of  frequency.  In 
hover,  for  exaaple,  large  low-frequency  inputs  are 
not  possible.  High-frequency  input  amplitudes  are 
also  often  reduced  because  rotor  or  structural 
resonances  and  actuator  authority  Units  are 
encountered.  Thus,  nonlinear  characteristics  nay 
not  be  exposed  consistently  across  the  frequency 
range,  or  froa  one  repeat  run  to  the  next.  In 
general,  these  inconsistencies  are  considered 
acceptable  since  the  pilot's  input  aaplitude 
during  the  frequency  sweeps  is  representative  of 
his  noraal  operating  control  technique.  The 
resulting  nodels  are  describing-functions  which 
reflect  the  "average"  dynamics  of  the  rotororaft 
in  its  noraal  operating  envelope . 

While  the  sweep  can  be  caaputer-generated  (as 
in  Ref.  40),  pilot-generated  inputs  have  aaae 
inportant  advantages: 

(1)  Since  pilot  test  technique  differs  aoae- 
what  in  each  repeat  run,  the  concatenated  records 
generally  contain  good  spectral  content  over  the 
entire  frequency  range  of  Interest.  The  differing 
input  amplitudes  for  each  run  yield  a  aodel  which 
reflects  the  average  vehicle  dynamics  over  a  wide 
operating  range  of  amplitude  and  input  forms. 

(2)  Although  oouputer-ganerated  sweep  inputs 
are  symmetric,  the  responses  will  not  be  symmetric 
and  oan  cause  the  aircraft  to  drift  significantly 
away  from  the  reference  trim  oondition.  This 
behavior  is  apparent  in  the  computer-generated 
sweeps  of  Ref.  40.  When  the  pilot  generates  the 
sweep  input,  low-frequency  trim  inputs  oan  be  used 
to  correct  the  drifts  in  the  reference  condition. 

(3)  Finally,  the  manual  sweeps  give  the 
pilot  full  control  over  the  oonduet  of  the  tests 
and  the  magnitudes  of  the  inputs  which  is  a  par¬ 
ticular  concern  in  hovering  tests. 

Some  useful  guidelines  were  compiled  by  the 
pilots  who  executed  the  frequency-sweep  tests  on 
the  teetering  rotor  helicopter.™  They  are  given 
in  the  Appendix. 


Rotororaft-Parameter  Identification 

Dynamics  identification  methodologies  gener¬ 
ally  fall  into  two  categories:  frequency-domain 
and  time-domain.  Each  approach  has  its  Inherent 
strengths  and  weaknesses  which  make  it  best  suited 
for  particular  applications.  Time-domain  (maximum 
likelihood)  identification  (Fig.  17)  uses  Kalaan- 


filter  technology  to  obtain  unbiased  estimates  of 
state-space  model  parameters.  Such  a  aodel  facil¬ 
itates  a  direct  comparison  between  stability 
derivatives  obtained  in  the  wind  tunnel  and  those 
of  the  actual  flight  vehicle.  Transfer  functions 
and  frequency  responses  are  by-products  of  time- 
domain  identification  since  they  are  calculated 
from  the  extracted  state-space  model  but  are  not 
identified  explicitly.  Thus,  nodels  which  provide 
a  good  fit  in  the  tine-domain  do  not  necessarily 
yield  transfer  functions  that  accurately  reflect 
the  frequency-domain  behavior.  One  reason  for 
this  is  that  time-domain  identification  uses 
least-squares  fitting  along  a  linear  time  scale, 
while  transfer  functions  and  frequency  responses 
are  displayed  on  a  logarithmic  scale  of  fre¬ 
quency.  Therefore,  time-domain  identification 
techniques  weight  their  results  more  heavily  at 
low  frequency  where  the  majority  of  the  data 
points  are  concentrated. 

While  time-domain  identification  methods  have 
been  extensively  applied  to  fixed-wing  vehicles, 
similar  experience  with  rotororaft  is  limited.  In 
the  Federal  Republic  of  Germany,  the  German  Aero¬ 
space  Research  Establishment  (DFVLR)  has  conducted 
extensive  studies  in  parameter  Identification  of 
rotororaft  dynamics  using  maximum  likelihood  time- 
domain  techniques.  Much  of  this  research  has  been 
associated  with  the  highly  coupled  BO- 105  hinge¬ 
less  rotor  helicopter. 20,21  In  the  national 
Research  Council  of  Canada,  a  time-domain,  quasl- 
1 inear  response-error  method  was  used  for  identi¬ 
fication  of  helicopter  stability  and  control 
derivatives.211  Recent  work  in  time-domain  identi¬ 
fication  of  rotororaft  has  also  been  reported  by 
researchers  at  the  Royal  Aircraft  Establishment 
(RAE),  Bedford,  England. 2? 

The  frequency-domain-identification  approach 
depicted  in  Fig.  18  uses  spectral  methods  to 
extraot  the  frequency  responses  between  selected 
input  and  output  pairs.  The  identification 
results  are  presented  in  Bode-plot  format:  magni¬ 
tude  and  phase  versus  frequency.  These  nonpara- 
metric  identification  results  are  very  useful  for 
flight-control  system  design  and  handling  quali¬ 
ties  studies.  Currently  proposed  handling- 
qualities  criteria  for  the  LHX  (Ref.  41)  are  based 
on  frequency-domain  parameters  which  can  be  read 
direotly  from  these  graphical  results.  Frequency 
responses  obtained  from  real  time  and  nonreal  time 
simulations  can  be  compared  directly  with  the 
flight  data  to  expose  limitations  and  discrepan¬ 
cies  in  the  simulator  models.  The  fact  that  this 
comparison  can  be  made  without  an  a  priori  assump¬ 
tion  of  model  structure  or  order  is  especially 
important  for  verifying  mathematical  models  of  new 
airoraft  configurations.  Tabulated  frequency 
responses  are  fitted  with  analytical  transfer- 
function  models  to  extraot  modal  characteristics 


which  are  useful  for  handling-quality  specifica¬ 
tions  given  in  lower-order  equivalent  system  terms 
and  for  transfer  function-based  control  system 
design  studies.  Since  this  fitting  procedure  is 
completed  after  the  frequency  response  is 
extracted,  the  order  of  the  transfer  function  can 
be  carefully  selected  to  avoid  an  overparameter- 
ixed  model.  Multi-input/multi-output  frequency- 
response  methods  are  suitable  for  extracting  a 
transfer  matrix  which  includes  the  Important  cou¬ 
pling  effects. 3° 

Once  the  transfer  functions  are  identified, 
state-space  models  given  in  terms  of  modal  or 
canonical  coordinates  can  be  realized.  State- 
space  models  can  also  be  extracted  in  terms  of  the 
physical  state  variables  using  least  squares  fit¬ 
ting  of  the  transfer-function  parameters  or  fre¬ 
quency  responses. 28  Finally,  the  extracted  models 
are  driven  with  the  flight  data  to  verify  the 
time-domain  characteristics.  The  semilog  fre¬ 
quency  format  of  the  Bode-plot  presentation  and 
subsequent  transfer-function  fit  makes  the  identi¬ 
fied  transfer-function  and  state-space  models  most 
accurate  at  mid  and  high  frequency  (initial  time 
history  transients).  The  low-frequency  and 
steady-state  response  prediction  of  the  extracted 
models  is  generally  not  as  good  as  in  the  time- 
domain  identification  approach. 

Miile  frequency-domain  identification  of 
fixed-wing  aircraft  has  been  conducted  for  many 
years,  its  application  to  rotorcraft  identifica¬ 
tion  is  fairly  recent: 

1.  XV-15  tilt  rotor  aircraft.29,30 

2.  Bell -214-ST. 39 

3.  HASA/Army  CH-47B.18 

4.  Canadian  MRC-205.40 

The  first  three  aircraft  are  being  investigated  by 
the  authors  and  are  discussed  in  the  following 
sections. 

Frequency-Domain  Identification  of  the  XV-15 

The  frequency-domain  identification  of  the 
XV-15  tilt-rotor  aircraft  (Fig.  19)  is  an  exten¬ 
sive  on-going  research  effort.  The  objectives  of 
this  project,  begun  in  1983,  are  to  identify  and 
document  the  open-loop  dynamics  of  the  XV-15 
throughout  its  flight  envelope,  to  compare  the 
simulation  and  flight  responses,  and  to  develop 
and  verify  a  transfer-function  model  description 
for  future  use  in  control -system  development. 
Identification  methodologies,  computational  tech¬ 
niques,  and  discussion  of  results  are  covered  in 
detail  in  Refs.  29  and  30.  Included  in  this 
research  have  been  multi-input/multi-output 


identification,  multi-axis  transfer-function 
fitting,  and  structural-dynamics  identification. 
The  current  work  involves  the  comparison  of 
frequency-domain  and  time-domain  identification 
results  for  the  hover  condition. 

Identification  of  Lower-Order  Equivalent  Transfer 
Functions  of  Bell  214-ST 

Frequency-domain  testing  of  the  Bell  214-ST 
single-rotor  helicopter  (Fig.  20)  was  completed  in 
October  1985  in  support  of  the  Army's  development 
of  an  updated  Mil-H-8501  and  LHX  handling- 
qualities  specification.  '  These  tests  were 
intended  to  demonstrate  the  frequency-domain  iden¬ 
tification  procedure  on  a  single-rotor  helicopter 
since  much  of  the  previous  methodology  development 
had  been  completed  on  the  (twin-rotor)  XV-15  air¬ 
craft.  A  discussion  of  the  testing  procedures  and 
Identification  results  for  the  Bell-214-ST  is 
presented  in  Ref.  39.  This  project  showed  the 
feasibility  of  frequency-sweep  testing  and  fre¬ 
quency-domain  identification  for  use  in  helicopter 
handling-qualities  specifications.  Further,  the 
results  showed  that  very  low  order  transfer- 
function  models  can  accurately  predict  the  large 
motion  time-domain  behavior  of  helicopters— even 
in  hover.  An  example  of  the  results  excerpted 
from  the  report  on  the  Bell-214-ST  test39  is  given 
below  to  illustrate  the  frequency-domain  identifi¬ 
cation  approach. 

The  lateral  stick  frequency-sweep  input  and 
the  corresponding  input  auto  spectrum  for  the 
hover-flight  condition  (with  the  control  system 
engaged)  was  previously  presented  in  Figs.  15 
and  16.  The  associated  roll  rate  response  for 
these  two  frequency  sweeps  is  shown  in  Fig.  21. 

The  maximum  roll  rate  is  about  ±15°  per  sec,  with 
somewhat  lower  values  for  very  low  and  very  high 
frequency  inputs.  The  corresponding  output  auto 
spectrum  (Fig.  22)  shows  that  the  roll  rate  exci¬ 
tation  is  roughly  constant  in  the  frequency  range 
of  0.3  to  7.0  rad/sec,  which  corresponds  to  the 
frequency  range  over  which  the  input  auto  spectrum 
is  constant  (Fig.  16).  At  high  frequency,  the 
spectral  response  drops  off  owing  to  the  K/s 
response  roll-off  of  the  rigid  body  helicopter 
dynamics  and  to  the  reduced  pilot  inputs  at  high 
frequencies  (Fig.  15).  The  peak  which  occurs  at 
u  =  11.9  rad/sec  is  due  to  the  excitation  of  the 
rotor-pylon  structural  mode.  The  output  auto 
spectrum  drops  sharply  for  frequencies  below 
0.1  rad/sec  since  there  is  very  little  pilot  input 
at  these  low  frequencies. 

The  frequency  response  of  roll  rate  to  lat¬ 
eral  stick  (P/6{,at)  la  shown  in  Fig.  23.  At  the 
higher  frequencies  considered,  the  response  exhib¬ 
its  a  K/s  characteristic  which  Indicates  that  a 
constant  roll  acceleration  results  from  an  Initial 


(4) 


input  of  lateral  stick.  The  presence  of  the 
rotor-pylon  mode  at  u  =  11.9  rad/sec  is  clearly 
seen  in  the  figure.  At  very  low  frequencies,  the 
roll  rate  is  reduced  significantly  because  of  the 
large  lateral  velocity  perturbations. 

The  quality  of  the  identified  frequency 
response  can  be  assessed  from  the  coherence 
function  shown  in  Fig.  24.  This  frequency- 
dependent  parameter  may  be  interpreted  as  that 
fraction  of  the  output  spectrum  estimate  which  can 
be  accounted  for  by  linear  relation  with  the  input 
spectrum  estimate.  Thus,  when  the  process  under 
investigation  is  perfectly  linear  and  the  spectral 
estimates  are  noise  free,  the  coherence  function 
will  be  unity  for  all  frequencies  in  the  excited 
input  spectrum  range.  A  value  of  the  coherence 
function  less  than  unity  will  result  from 
nonlinearities  in  the  system,  input/output  noise, 
or  cross-coupled  control  inputs.  As  can  be  noted 
in  Fig.  24,  good  frequency-response  identification 
is  achieved  in  the  range  of  0.2  to  12.0  rad/sec. 
The  oscillation  in  the  coherence  function  for 
frequencies  greater  than  12.0  rad/sec  and  the 
rapid  decline  in  coherence  function  for  fre¬ 
quencies  below  0.3  rad/sec  are  strong  indications 
of  reduced  spectral  estimate  accuracy.  The 
identification  at  low  frequency  can  be  improved  by 
reducing  the  period  of  the  first  two  frequency 
sweep  cycles  (Fig.  15).  High-frequency  excitation 
for  the  Bell-214-ST  tests  was  restricted  because 
of  the  excitation  of  the  rotor  pylon  mode. 

The  next  step  in  the  identification  procedure 
is  the  fitting  of  lower-order  transfer-function 
models.  In  the  case  of  the  XV-15  study,  the  open- 
loop  coupled  dynamics  were  sought.  Therefore,  the 
selected  model  orders  corresponded  to  the  physical 
order  of  the  coupled  open-loop  dynamics. 3°  Thi3 
is  consistent  with  the  earlier  discussion  on  model 
order  selection  in  the  present  paper.  In  the 
Bell-214-ST  study,  closed-loop  transfer-function 
models  were  sought.  The  overall  effect  of  the 
closed-loop  augmentation  is  to  drive  the  open-loop 
poles  into  the  zeros,  thereby  suppressing  the 
open-loop  inter-axis  coupling  and  reducing  the 
residues  of  some  of  the  open-loop  modes.  An 
appropriate  clo3ed-loop  representation  for  the 
purposes  of  thi3  study  was  one  that  accurately 
modeled  the  dominant  on-axis  closed-loop  responses 
and  ignored  the  off-axis  responses  and  nearly 
canceled  modes.  Therefore,  the  transfer-function 
order  was  selected  as  a  minimum  order  which  would 
yield  a  reasonable  fit  of  the  flight  data.  In  the 
case  of  the  roll  (and  pitch)  closed-loop  attitude 
dynamics  of  the  Bell-214-ST,  the  following  second- 
order  transfer-function  form  was  found  to  ade¬ 
quately  reflect  the  dynamics  of  the  helicopter 
both  in  hover  and  cruise: 


(s) 


Ks  e 


-IS 


LAT 


(3  +  l/T^fs  +  1/T2) 


The  exponential  term  in  the  numerator  is  included 
to  account  for  the  phase  lag  due  to  high-frequency 
dynamics  such  as  the  rotor  and  actuators. 

The  maximum  frequency  range  usable  for 
transfer-function  fitting  corresponds  to  the 
region  of  maximum  coherence  (0.2-12.0  rad/sec)  in 
Fig.  24.  However,  the  2nd-order  model  of  Eq.  (4) 
only  accounts  for  the  rigid  body  dynamics,  so  the 
fit  must  be  restricted  to  the  frequency  range 
which  excludes  the  rotor-pylon  resonance.  Thus, 
the  "frequency  range  of  model  applicability"  is 
0.2  to  10.0  rad/sec.  This  approach  is  consistent 
with  the  earlier  discussion  of  vertical 
acceleration  models. 


Applying  the  transfer-function  form  of 
Eq.  (4)  to  fit  the  roll  rate  response  of  Fig.  23 
in  the  range  of  0.2  to  10.0  rad/sec  yields  the 
following  numerical  parameters: 


K  s  35.41  deg/sec2/ in.  ^ 
x  s  0.050  sec 


i-  =  0.38  rad /sec 
l1 


>  (5) 


2.83  rad/sec  J 

i  Eq.  (5),  it  can  readily  be  seen 
that  the  roll  response  is  heavily  damped  (no 
complex-valued  roots)  with  a  dominant  response 
time  constant  of  T?  =  0.3  sec.  The  small  effec¬ 
tive  time  delay  of  x  =  0.050  sec  shows  that  the 
rotor,  actuator,  and  structural  dynamics  do  not 
cause  significant  phase  lag  in  the  frequency  rate 
of  concern  for  piloted  handling  qualities  (about 
0.1  to  10.0  rad/sec).  As  shown  in  Fig.  23,  the 
extracted  model  fits  the  data  very  well  in  the 
frequency  range  of  applicability  (0.2-10  rad/sec). 


T_ 


Referring  to 


The  ability  of  this  simple  transfer-function 
model  to  predict  the  time-domain  roll  response 
characteristics  in  hover  was  tested  using  aircraft 
response  to  step  inputs.  This  input  was  selected 
for  the  verification  study  to  show  the  robustness 
of  the  extracted  model  in  predicting  response 
characteristics  to  input  forms  other  than  the 
frequency-sweep  form  used  in  the  identification 
process.  However,  as  discussed  earlier,  the  step 
input  has  high-frequency  spectral  content  which 
excites  aircraft  modes  outside  of  the  transfer- 
function  model's  range  of  applicability 
(0.2-10.0  rad/sec).  Most  Important  is  the 
unwanted  excitation  of  the  rotor-pylon  resonance 


(«  *  11.9  rad/aec).  Therefore,  the  Input  and 
output  flight  data  are  low-pass  filtered  to  atten¬ 
uate  the  speetral  content  for  frequencies  beyond 
«  »  10.0  rad/sec.  The  filtered  step  Input  flight 
data  are  used  to  drive  the  transfer-function  model 
for  comparison  with  the  filtered  output  flight 
data.  These  filtered  time  histories  are  shown  In 
Fig.  25. 

A  fairly  large  step  input  of  «LAT  =  0.75  in. 
was  applied  and  held  constant  until  a  maximum 
steady-state  roll  rate  was  achieved  (about  7® /sec 
as  seen  in  Fig.  25).  The  transfer-function  model 
is  seen  to  closely  reflect  the  filtered  roll  rate 
response  over  the  entire  time  history  including 
the  recovery  phase.  The  differences  which  are 
apparent  in  the  regions  of  maximum  roll  rate 
mostly  result  from  using  a  single-axis  model 
(Eq.  (4)),  which  ignores  the  contribution  of  the 
off-axis  inputs.  The  very  good  overall  match 
between  the  flight  data  and  the  transfer-function 
model  response  shows  the  utility  of  this  fairly 
low-order  model  in  characterizing  the  closed-loop 
dynamics  of  the  Bell-214-ST.  The  fact  that  the 
prediction  is  very  good  for  this  (smoothed)  step 
input  form,  which  is  markedly  different  in  shape 
from  the  frequency  sweep,  shows  the  robustness  of 
this  identification  approach  and  the  advantage  of 
using  minimum-order  transfer-function  models. 

Identification  of  CH-47B  Vertical  Dynamics  Model 
In  Hover 

The  same  frequency-domain  procedure  used  in 
the  identification  of  the  two  preceding  airoraft 
was  also  applied  to  the  identification  of  the 
vertical  dynamics  of  the  CH-47B  research  aircraft 
(Fig.  26)  in  hover.  A  varying  frequency  input  in 
collective  pitch  was  employed  to  identify  first 
the  freouency  response  of  the  vertical  accelera¬ 
tion  to  collective  input.  Figure  27  shows  a 
sample  of  typical  frequency  sweeps  of  the  input 
with  frequency  varying  from  about  0.05  Hz  to  about 
2  Hz.  The  normal  acceleration  response  with  a 
discernible  lead  in  low  frequency  is  also  shown  in 
the  figure.  Because  of  some  sharp  edges  in  the 
frequency  sweeps,  substantial  power  in  the  collec¬ 
tive  input  beyond  2  Hz  (slightly  above  3  Hz)  is 
present,  as  can  be  seen  in  the  input  auto-spectrum 
and  input-output  cross-spectrum  plots  shown  in 
Fig.  26. 

The  extracted  frequency  response  as  shown  in 
Fig.  29  and  the  accompanying  coherency  function 
(Fig.  30)  Indicate  that  the  nonparametric  identi¬ 
fication  is  good  between  the  frequency  range  of 
0.2  to  20  r/s.  Hie  transfer-function  amplitude 
plot  shows  that  there  is  a  resonant  peak  at  around 
17  r/s  as  predicted  earlier  in  Fig.  7.  The  pre¬ 
dicted  phase  lead  in  the  low-frequency  region 


shown  in  Fig.  7  also  appears  in  the  identified 
transfer-function  phase  plot  in  Fig.  29. 

Guided  by  the  analytical  modeling  efforts 
discussed  in  the  previous  section,  a  series  of 
parametric  transfer-function  models,  ranging  from 
1st  to  5th  order  was  fit  to  the  identified  fre¬ 
quency-response  plots.  The  results  which  are 
listed  in  Table  7  are  classified  according  to  the 
order  of  the  denominator.  The  Ist-order  model 
corresponds  to  the  quasi-static  model  shown  in 
Table  3.  The  next  higher-order  model  selected  for 
fitting  the  frequency-response  plots  is  a  3rd- 
order  model  which  corresponds  to  the  no  dynamic- 
inflow  case  (Case  (b)  in  Table  3);  no  2nd-order 
models  were  used.  Two  4th-order  models  were 
used:  one  with  4th  order  in  numerator  (which 
corresponds  to  Case  (a)  in  Table  3)  and  the  other 
with  3rd  order  in  the  numerator  (which  corresponds 
to  the  case  having  no  direct  control  effective  on 
the  vertical  acceleration  due  to  blade  pitch 
changes) .  °  The  5th-order  model  was  used  to 
account  for  possible  effects  due  to  variation  in 
rotor  RPM. 

The  frequency-response  fits  for  a  1st-  and  a 
4th-order  model  (models  (a)  and  (e)  in  Table  7) 
are  shown  in  Fig.  31.  For  the  Ist-order  model,  it 
proved  to  be  difficult  to  fit  the  model  accurately 
over  the  entire  frequency  range;  the  fit  was, 
therefore,  limited  to  the  range  of  0.1  to  3  r/s, 
instead  of  0.1  to  20  r/s  as  applied  to  all  other 
models.  This  causes  the  Ist-order  model  to  match 
the  flight  data  better  at  low  frequency  than  the 
higher-order  models,  as  can  be  seen  in  Fig.  31. 

All  the  higher-order  models  exhibit  a  nonminimum 
phase  characteristic  as  predicted  in  the  analyti¬ 
cal  modeling  previously  shown  in  Table  3.  The 
index  for  the  fitting  error  improves  as  the  order 
of  the  transfer-function  numerator  and  denominator 
increases  because  of  more  free  parameters  avail¬ 
able  for  fitting  the  data.  It  is  noted  in 
Ref.  42,  however,  that  a  model  with  more  free 
parameters,  though  in  general  achieving  a  better 
fit  to  a  set  of  given  data,  does  not  necessarily 
provide  a  better  predictive  capability  (which, 
after  all,  is  the  very  purpose  of  developing  a 
model). 

As  a  means  for  testing  the  predictive  capa¬ 
bility  of  those  six  identified  models,  the 
vertical-acceleration  responses  to  a  step  input  in 
collective  input  (0.62  in.)  were  generated  from 
the  models.  (In  calculating  these  responses,  the 
small  time  delay  (or  advance)  contained  in  the 
models  was  neglected.)  These  responses  of  the 
identified  models  (Fig.  32)  compare  very  favorably 
with  those  of  the  analytical  model  (Fig.  6).  The 
responses  of  the  two  identified  Ist-order  models 
(a  and  b),  a  4th-order  model  (e),  and  the  5th- 
order  model  (f)  are  plotted  together  with  flight 


data  (of  Fig.  5)  In  Fig.  33.  As  expected,  the 
Ist-order  model  (a),  because  of  its  limited  fre¬ 
quency  range  of  applicability,  fails  completely  to 
predict  the  overshoot  of  the  vertical  acceleration 
owing  to  the  high-order  effects  of  inflow  and 
flapping  dynamics.  However,  the  long-term 
response  matches  well  with  the  flight  data.  The 
other  identified  Ist-order  model  (model  (b)), 
which  fits  over  the  entire  frequency  range  with  a 
large  fitting  error,  produces  a  response  that 
attempts  to  compromise  the  initial  and  final 
responses  of  the  flight  data.  As  a  result,  its 
response  is  inconsistent  with  the  responses  of 
other  models  and  its  range  of  applicability  and 
usefulness  becomes  rather  uncertain.  The  4th- 
order  model  provides  a  much  better  predictive 
capability  than  does  the  Ist-order  ones  for  the 
initial  transient.  The  response  of  the  5th-order 
model  matches  even  better  with  the  flight  data  in 
the  initial  transient,  but  its  long-term  response 
becomes  somewhat  worse  than  that  of  the  4th-order 
model,  indicating  that  overparameterization  may  be 
present  in  the  $th-order  model. 


Concluding  Remarks 

This  paper  has  discussed  in  some  detail  the 
Importance  of  recognizing  that  each  lower-order 
model  used  for  rotorcraft-parameter  identification 
has  a  limited  range  of  applicability.  Awareness 
of  this  fundamental  limitation  is  of  paramount 
Importance  in  designing  flight-test  experiments, 
in  processing  the  test  data,  and  in  interpreting 
the  identification  results.  Test-input  signals 
must  not  only  be  designed  to  optimize  the  informa¬ 
tion  content  in  the  rotorcraft  response  measure¬ 
ments  but  also  be  properly  conditioned  such  that 
the  highest  frequency  content  of  the  flight-test 
input*  is  consistent  with  the  frequency  range  of 
applicability  of  the  model.  Examples  have  also 
been  given  to  illustrate  the  potentials  of  using 
multi-axis  test  inputs  to  enhance  the  parameter 
identlfiabllity. 

Also  discussed  in  great  length  are  the  bene¬ 
fits  of  using  frequency  sweeps  as  input-test  sig¬ 
nals  for  identification  of  frequency-response  for 
the  rotorcraft  and  for  subsequent  fitting  of  para¬ 
metric  transfer-function  models.  The  frequency 
sweep  provides  a  good  control  over  the  frequency 
range  of  Interest  and  it  lends  Itself  naturally  to 
non parametric  frequency-domain  identification. 
Analytical  modeling  and  understanding  of  the 
physics  Involved  in  rotorcraft  flight  dynamics, 
especially  understanding  the  limitation  of  lower- 
order  models,  can  be  more  important  than  merely 
relying  on  the  identification  algorithms  in  the 
final  stage  of  fitting  parametric  transfer- 
function  models  to  the  flight-extracted  frequency 
response  data. 


Appendix:  Pilot  Comments  on  Frequency  Sweep 
Inputs  in  BHT1  214ST  Helicopter 

Pilots:  CW4  Robert  A.  Williams  and  Cpt.  Randy 
Cason,  USAAEFA  Edwards  AFB,  CA 

Two  Flights  Completed  2  August  1985 
Hover  Flight  1.2  hours 
Level  Forward  Flight  90  KCAS  1.6  hours 

Recommendations 

A.  Recommend  that  yaw  inputs  be  performed  prior 
to  other  inputs  since  the  yaw  rates  were  very 
low  and  it  will  aid  in  the  pilot  learning 
curve.  The  test  should  be  done  in  order  of 
increasing  difficulty:  yaw,  collective, 
longitudinal,  and  then  lateral  inputs. 

B.  Recommend  that  aircraft  gross  weight  be  at  a 
minimum  during  the  collective  inputs.  Maxi¬ 
mum  gross  weight  for  the  Bell  214ST  is 
17,500  lb  and  the  test  was  conducted  at 
approximately  13,000  lb.  Torque  readings 
were  consistently  greater  than  90%  (above 
100%  constitutes  an  overtorque). 

C.  Recommend  that  the  pilot  be  "coached"  during 
the  input.  It  is  very  easy  to  remain  at  one 
frequency  too  long,  having  the  engineer  tell 
the  pilot  to  dwell  on  a  specific  frequency 
longer  or  increase  frequency  during  a  data 
run  aided  data  acquisition.  This  assumes  the 
engineer  has  real  time  data. 

D.  Recommend  that  the  copilot  or  flight  test 
engineer  coach  the  pilot  for  the  low  fre¬ 
quency  responses  by  counting  seconds  for 
timing  the  quarter  period.  This  should  only 
be  done  for  the  lowest  frequencies.  It  was 
found  that  if  the  copilot  tried  counting  at 
the  higher  frequency,  it  only  mixed  up  the 
pilot  and  resulted  in  the  pilot  following  the 
copilot's  counting  rather  than  increasing  the 
frequency  as  required  by  the  test. 

E.  Recommend  that  pilots  practice  the  inputs 
utilizing  external  power  supply  (hydraulic 
and  electrical  power)  on  the  ground  prior  to 
testing.  This  should  significantly  reduce 
test  flight  time. 

F.  Recommend  a  "scope"  with  "freeze"  capability 
be  installed  in  the  test  aircraft  unless 
"real  time"  TM  is  available. 

G.  Recommend  that  the  test  team  be  briefed  on 
all  "aircraft"  natural  frequencies  below 
4.0  Hz  and  any  resultant  problems  which  may 
be  encountered  at  these  frequencies. 


H.  During  flights  that  require  longitudinal 
inputs  the  Instruaented  A/S  booa  oust  be 
monitored  olosely  for  deflection  beyond 
limits  or  removed. 
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Table  la  F  and  G  aatrloea  for  10th  order  teetering  rotor  hel loop ter  la  hover 
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-0.001641 

13.14 

21.65 
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State  vector  *  (a1  b,  a1  b,  •  *  q  p  u  v)' 
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Table  1b  F  and  C  eat rices  for  simplified,  quasi-static 
6th  order  helloopter  Model 
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Table  3  Vertical  acceleration  to  collective  transfer-functions  of  the  three  vertioal 

dynaalc  aodels 


a.  *»  .  -3  9*  -  22  a3  ♦  8771  a2  ♦  47630  s _ 

*c  s4  ♦  34.3  a3  ♦  698.3  a2  ♦  5556  a  ♦  1520.3 
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Table  4  Comparison  of  allaroo  Inputs  and  ruddar  inputs 
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Table  5  Suboptiaal  inputs  -  simultaneous  use  of  aileron  and  rudder 
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Table  6  Reeoawnded  practical  aileron  and  rudder  Input 
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Fi* .  1  Comparison  of  tho  frequency  responses  of  P/*1  for  10th  and  6th  order  sysi 
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Fig.  2  Coopori— a  of  tho  fraquancy 


Pig.  4  Comparison  of  tho  transient  responses  of  the  10th  and  6th  order  Models  to  a  high 

frequency  (9  rad/sec)  doublet  input. 


Fig.  5  Vertical  acceleration  response  of  the  CH-47B  helicopter  to  a  step  collective  Input, 

i«c  z  0.62  in.  (A»q  z  0.0201  rad). 
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Fig.  6  Calculated  vertical  acceleration  response  to  0.62  in.  (0.0201  rad)  step  collective 
input:  (a)  4th -order  nodel,  (b)  3rd -order  model,  (c)  Ist-order,  quasi-static  model. 
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Pig.  8(b)  Tim  history  and  frequency  distribution  of  s  aaooth  test  input  end  s  nultistep 

test  input. 
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Pig.  9  Sohemtle  diagram  for  Input  signs!  conditioning 
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Fig.  10  Cooparlson  of  the  transient  responses  of  the  10th  and  6th  order  Models  to  a  high 
frequenoy  (9  rad/sec)  doublet  shaped  with  a  2nd  order  filter  (wc  «  2.5  rad/see,  c  >  0.7). 
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Fig.  It  Aircraft  response  to  a  rudder  doublet  input. 


rtf.  22  Moll  rate  output  autospectrua  (Boll  214-ST). 


rig.  23  Boll  attitude  roopouao  to  lateral  stick  In  hover  (Boll  214-ST). 


flf.  24  Coherence  function  for  roll  rate  response  identification  (Bell  214-ST) 
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IS 


Fig.  25  Comparison  of  filtered  aircraft  response  and  transfer  function  model  response  to 
filtered-step  input  (Bell  2H-ST):  (a)  lateral  stick  deflection,  (b)  roll-rate  response. 


Fig.  26  CH-47B  variable-stability  research  helicopter. 


Pig,  28  Auto>apeetrua  of  collective  Input  and  croaa-apectrun  of  collective  input  and 
vertioal  acceleration  reaponae  (CH-47B  flight  data). 


FREQUENCY,  rad/we 

Pig.  29  Frequency  responses  of  vertical  acceleration  to  collective  input  (CH-47B  flight 


FREQUENCY,  rad/aoc 

Coherence  function  associated  with  transfer-function  identification. 
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Fig.  31  Transfer- function  aodels  fit  to  flight  extracted  frequency-response  data 

(a)  aagnitude,  (b)  phase. 
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Fig.  32  Predicted  ^  responses  to  step  collective  input  (0.62  In.)  from  the  six 

identified  models. 


Pig.  33  Flight  date,  predicted  responses  of  1st,  4th,  and  5th  order  models. 


